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ABSTRACT 


Three-dimensional trusses are designed for minimum 
weight, subject to constraints on: member stresses, Euler 
buckling, joint displacements and system natural frequencies. 
Multiple static load conditions are considered. 

The finite element displacement method of analysis is 
used and eigenvalues are calculated using the subspace 
iteration technique. All gradient information 1S calculated 
analytically. 

The design problem is cast as a multi-level numerical 
Optimization problem. The joint coordinates are treated as 
system variables. For each proposed configuration, the 
member sizes are updated as a sub-optimization problem. 

This sub-problem is efficiently solved uSing approximation 
concepts in the reciprocal variable space. The multi-level 
approach is shown to be an effective technique which con- 
veniently takes advantage of the most efficient methods 
available for the member sizing problem. 

Examples are presented to demonstrate the method. The 
optimum configuration is shown to be strongly dependent on 


the constraints which are imposed on the design. 





TAB EE OR eon ehilsS 


a. INTRODUCT ION--------------- ------------------------- eZ 
et MATHEMATICAL FORMULATION---------------------------- 16 
A. INTRODUCTION------- 3-3-2 2 ee 16 

B. ANALYSIS 7-339 ---  - 1G 

1. Static Analysis----------------------------- ak 

2. Dynamic AnalysiS------------------ = ie: 

C. ANALYTIC GRADIENTS OF THE CONSTRAINTS ----------- 1S 


1. Gradient of Member Stresses with Respect to 
the Reciprocal of Area Variable------------- 20) 


2. Gradient of Nodal Displacements with Respect 
to the Reciprocal of Area Variable---------- 2 


3. Gradient of Frequencies with Respect to the 
Reciprocal of Area Variable----------------- 22 


4, Gradient of Stresses with Respect to the 
Joint Coordinate VariableS------------------ 23 


5. Gradient of Displacements with Respect to 
Powe eCOchamMictem Vania les ————————————— = oe 


6. Gradient of the Natural Frequencies with 


Respect to the Joint Coordinate Variables---24 
APPROXIMATION CONCEPTS-------------------------- 25 
OBJECTIVE FUNCTION------------------------------ 27 
CONSTRAINTS----------- ~------------------------- 27 
1 BeOS S=— eS = ee —- oa Zi 
Pe nee ee abn = Se 28 
Sess clCy)—-————aeaee = =e oe So SSeS 28 





ee ee SO nena aS = = — — — — a 28 


5. Displacement-------------------------------- 28 

6. Limits on Coordinate Variables-------------- 29 

G. GENERAL FORMULATION----------------------------- Zo 

Ill. OPTIMIZATION---------------------------------------- og 
A. INTRODUCTION+- == -- 7 - - = - - - - - - - - - - - - - - - - - - -- 31: 

B. GENERAL FORMULATION-<---------------------------- ort 

C. OPTIMUM GEOMETRY DES IGN--- 9-7 9-9 ern - BV 

Pee we GuoMarn:, DES TGN==—=—<———-————-———==—-——————— 34 

iy. NUMERICAL EXAMPLES rr rrr nn nr nr nr rr rrr rrr rr ree 36 
A. INTRODUCTION or or nr nr rn rr rrr rrr rrr eee 36 

B. CASE 1: 18-BAR PLANAR TRUSS-------------------- 36 

Ll. CaS@ Larner rrr rr rr rr rrr re S7 

2. Case lb--- een en ne rere 38 

3. CaSe Lerten ener een nn rn rr rr reer 38 

ES oO ee ee 38 

C. CASE 2: 47-BAR PLANAR TOWER-------------------- 39 

D. CASE 3: 25-BAR SPACE TOWER--------------------- 40 

ie CONCLUSIONS AND RECOMMENDATIONS--------------------- 42 
A. CONCLUSIONS --------------------- 99-99 nr 42 

B. RECOMMENDATIONS--------------------------------- 43 
APPENDIX An--- 7-992 nn rr rrr rere 44 
LIST OF REFERENCES ------------ 9-9 enn nnn 70 
PiTCAn DISTRIBUTION LIST------------------------~--------- 72 








ive. 


oo 


i. 


Wales 


ert. 


weed. 


x. 


ral, 


er. 


Pert. 


XIV. 


LIST 


18-Bar Planar Truss, 


OF TABLES 


Loads and ConstantS------------ 50 


18-Bar Planar Truss, Design Information. Stress----- oa 
18-Bar Planar Truss, Design Information. Stress, 

Buler Buckling-------------------------------------- 58 
18-Bar Planar Truss, Design Information. Stress, 

Euler Buckling, Displacement------------------------ a2 
47-Bar Planar Tower. Load Conditions--------------- 60 
47-Bar Planar Tower. Load Conditions--------------- 61 
47-Bar Planar Tower. Area. Stress, Euler Buckling-62 
47-Bar Planar Tower Coordinates. Stress, Euler 
Ce = oo 
47-Bar Planar Tower. Coordinates. Stress, Euler 

Bere amen sO aCe Gen eS 64 
47-Bar Planar Tower. Areas. Stress, Euler Buckling, 
Pesce meh aw ent a Se a a 8 8 8 === s5 65 
47-Bar Planar Tower. Areas. Stress, Euler Buckling, 
Prcfemccenete, dua. frequency —————— == ———<<<<<=-<-=<--=- 66 
47-Bar Planar Tower. Coordinates. Stress, Euler 
BuGi@king, Displacement, and Frequency —---<<<<-<---=— 67 
Pe=pomecpace: Tower. Load Conditions—----==-<-------- 68 
Zo=pan space Tower. Design Inftormatlon------------- 69 





Eiol OF FLGURES 


Pee l-. Wop bawsoboanare cuss. oeress, Hudier Buckling, 
Displacement, and Frequency Constraint------------ 44 

eee 18-Bar Planar Truss. Weight vs. Iteration Number, 
StreS $299 - - - + + 45 

A 18-Bar Planar Truss. Weight vs. Iteration Number. 
Stress, Euler Buckling---------------------------- 46 


pee 4. 18-Bar Planar Truss. Weight vs. Iteration Number, 
Stress, Euler Buckling, Displacement-------------- 47 


Po 18-Bar Planar Truss. Weight vs. Iteration Number. 
Stress, Euler Buckling, Displacement, and 


Frequency ~~ - 7-9 - - - - - - - - - - - - - - - - - - - - - 48 
AO: 47-Bar Planar Tower---~--~--~------~--~--------~---—-—-—- 49 
Bees 47-Bar Planar Tower. Stress, Euler Buckling, 

Displacement, Frequency ~-------- 9-9 ee ee ee 50 


A.8. 47-Bar Planar Tower. Weight vs. Iteration Number. 


Stress, Euler Buckling-----37~----93-<23 3-2 ---------- a 
PES 47-Bar Planar Tower. Weight vs. Iteration Number. 
Stress, Euler Buckling, Displacement---~----------- oie 


A.10. 47-Bar Planar Tower. Weight vs. Iteration Number. 
Stress, Euler Buckling, Displacement, Frequency--- 53 


A.ll. 25-Bar Space Tower------------- 9999-2 = 54 


A.12. 25-Bar Space Tower. Weight vs. Iteration Number. 
Stress, Euler Buckling, Displacement, Frequency--- 25 





is is Cha ORR OHA 
Q ase 


= ws 
tt 


AC 


Za 


eC 


Oi Ni 


(A 


List OF SYMBOLS 


Cross-sectional area of member i 


Minimum and maximum allowable area variables, 
respectively 


Constants defined by Eq. 43 

Factordefined by Eq. 31 

Matrix of direction cosines 

Young's Modulus for member i 

Vector of applied loads 

Constraint function of coordinate variable 
Constraint function of area variable 
Active constraints of area variable 
Identity matrix 

Global stiffness matrix 

Element stiffness matrix 

Length of a bar 1 

Global mass matrix 

Global generalized mass matrix 

Lumped mass matrix 

Number of active constraints 

Number of initial constraints 

Direction vector 

Direction vector on the coordinate design space 


Direction vector on the area design space 





4 


a, U, V 


Coe, 2 


Nodal displacement vector 

ith nodal displacement vector 

Minimum and maximum allowable displacement value 
Vector of design variables 

Vector of coordinate design variables 

Vector of area design variables 


lth independent design vector 


Minimum and maximum allowable coordinate variable 
Dmeect Lomvecosines 


Scalar parameter defining distance of move in 
the area or coordinate design space. 


Objective function of direction function problem 
Displacement constraint value 

Lower and upper limits on displacement constraints 
Kronecker delta 

Push-off factor 

Stress in member 1 

Allowable compressive stress in member 1 

Allowable tensile stress in member i 

Allowable Euler buckling stress in member i 
Material density 

EBigenvector corresponding to Ws 

ith eigenvalue 

Minimum and maximum allowable fundamental frequency 


Nodal matrix 
Expectral matrix 
Frequency limit 


10 





ACKNOWLEDGEMENT 


The author wishes to express his deep appreciation to 
his faculty advisor, Professor Garret N. Vanderplaats, for 
his continuous guidance, assistance and encouragement as 
well as his great enthusiasm during the course of this 
project; to Professor David Salinas who read this project, 
and gave of his time and engineering expertise; to the 
distinguished members of the faculty and the staff of the 
Department of Mechanical Engineering for their support in 
very way, my Sincere appreciation. 

This investigation would not have been possible without 
mgemracrlities and support of the W. R. Church Computer 
C=RiEcr. 

Infinite thanks and grateful appreciation go to my wife, 
Yolanda, my son, Paulo, and daughter, Maria Gabriela, whose 
understanding and encouragement made the difficult time 


bearable. 


Ana 








Pen LRObUeTICN 


The process of optimization of structures has undergone 
important changes since its development in the early 1960s. 
Minimum weight of elastic truss structures subject to multi- 
ple loadings has been an active area of research. Attention 
has been focused on the problem of least weight, when the 
overall layout was known in advance, and when the cross- 
sectional area was the design variable. Some attention has 
been directed toward the optimum configuration of the 
Structure. DeSign improvements in this area often exceed 
those in fixed-geometry and so shape optimization is of 
major interest. 

Pioneer work in shape optimization was conducted in 1964 
[Ref. 1] by Dorn, Gomory and Greenberg. The optimal con- 
nectivity of nodes for truss members, subject to a single 
load condition, was found and minimum weight designs were 
achieved. In their work only planar trusses were tested, 
and the process was presented as a plastic design problem 
using linear programming. 

Their work was followed by Dobbs and Felton [Ref. 2] 
who in 1969 investigated the effect of multiple load condi- 
tions on the optimum configuration of trusses through the 


use of non-linear programming methods. Again, only planar 


HZ 








trusses were considered, subject to failure by stress and 
elastic buckling. 

Later, Pedersen [Ref. 3] considered the positions of 
eae jOints as continuous design variables in addition to 
the areas of the bars. Stress, displacement, and buckling 
constraints were considered. Pedersen's work is significant 
because the optimization process is carried out by consider- 
ing two separate design spaces. The optimization is 
achieved by successive iterations using a gradient method 
with move limits. 

The optimization process was advanced by Vanderplaats 
and Moses [Ref. 4] who divided the design space into two sub- 
Spaces, separating the area variables and the joint position 
variables. Multiple load conditions and constraints on 
stresses and Euler buckling were considered. The optimiza- 
tion was carried out alternatively between the two spaces 
until convergence was achieved. Three-dimensional indeter- 
Minate trusses were designed subject to multiple loading 
conditions. This work was extended by Vanderplaats [Ref. 5] 
to include displacement constraints. 

In other research, Spillers [Ref. 6] considered statically 
indeterminate trusses. The optimization followed an itera- 
tive design where the member sizes and node locations were 


the design variables. 


ds 





Recently, Imai [Ref. 7] treated the sizing and con- 
figuration variables simultaneously for either determinate 
or indeterminate trusses. The optimization was achieved 
using the Augmented Lagrange multiplier method. 

The design problem considered in this study is the opti- 

mum configuration of three-dimensional indeterminate 
trusses, for multiple prescribed static load conditions. 
The objective is to minimize the weight of the structure, 
design variables are the node coordinates and the member 
Sizes. Constraints inciude stress, Euler buckling, dis- 
placement, and the natural frequency of the system. 

Some approximation concepts are introduced in order to 
reduce computational effort and to reduce the nonlinear 
characteristic of some constraints. Among these, a first 
order Taylor Series expansion is applied to approximate the 
Constraints. 

The optimization proceeds iteratively in two design 
Spaces: the member sizing space, where the structure 1s 
optimized for a fix layout, and the coordinate space, where 
the geometry is allowed to vary. In both optimization 
processes the minimum weight of the structure was maintained, 
subject to the requirements that the constraints remain 
satisfied. 

The mathematical formulation is presented in Chapter II. 


The objective function and constraints are defined in terms 
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of the design variables. The analytic gradients are also 
formulated. 

The optimization technique is discussed in Chapter III. 
Several examples are presented and explained in Chapter IV. 
Conclusions and recommendations for future work are 


presented. 


ILS 








If. MATHEMATICAL FORMULATION 


or LN RODUCTION 

several features are desirable when finite element 
methods of structural analysis are used in optimization. 
First, the number of analyses for the structure should be 
kept as low as possible. Second,the amount of gradient in- 
formation required during the design process should be 


reduced as much as possible. 


Bee ANALYSIS 
1. Static Analysis 

The initial layout of the truss, the member sizing 
and material properties (which may be different for each 
member), a set of external loads, and support conditions 
are initially specified. 

The analysis for the stresses and deflections must 
be carried out satisfying the conditions of equilibrium of 
forces at the nodes and compatibility of deformation. If 
the material of the structure behaves in a linear manner, 
Hooke's law will establish the force-deflection relationship. 

For a truss, it 1s also necessary to establish the 
following assumptions before selecting the method Be comput- 
ing the internal forces. The discrete element is treated 


as pin-connected, and loads and reactions are supported at 
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the joints. In this study, the weight of the members is 
not included as loads. 
The Displacement (Stiffness) method [Ref. 8] 
considers the joint displacement components as the unknowns, 
and 1s written in the most general form using matrix notation, 


» 


su =e (Ga. =) 


roy 


where K is the global stiffness matrix, F is the vector or 
vectors of applied loads, and u is the vector or vectors of 
displacements. Equation 1 is the set of equilibrium equa- 
tions, and is formulated such that the compatibility is auto- 
matically satisfied. 

The generality of the method 1S important if either 
statically determinate or indeterminate trusses are analyzed. 
The global stiffness matrix is symmetric and sparce. These 
features are used to write the code for a computer solution, 
and the matrix K is stored in compact form for efficient 
numerical solution. 

Once the displacements at every node are known, the 


internal forces or stresses are calculated by applying force- 





deflection relations. This is defined as: 
Wie a Cee 
E. p? SG Kg | 
= a — . Eq. 2 
qj 5 ime Yo4 7 Yr3 ee 
yn kg 


where i and j are element and load condition numbers respec- 
tively, and k and 2 are node numbers associated with the 
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element 1. oa 1s the stress, Ee 1s the Young's modulus, 
and Le is the length. Matrix D = f(a,u,v) contains the 
direction cosines. For brevity, hereafter the second sub- 
Script 1s omitted and it is assumed that the stress or 
displacement corresponds to the appropriate loading condition. 
2. Dynamic Analysis 

When system natural frequency constraints are imposed 
in the design process, the corresponding dynamic analysis of 
the structure has to be carried out. This requires the 
solution of an eigen-problem to determine the natural fre- 
quencies and normal modes. For linear elastic structures, 
the finite element approach leads to the well-known equation 
SeemocLon, COnSidering free vibration SCOnclerons, 


Mu+ Ku = 0 (cea) 


where M is the global mass matrix, and u is the linear accel- 
eration. Assuming a solution of the form 

y= pe tet (Eq. 4) 
where w is the angular natural frequency of vibration of 
the structure, and 6 is the corresponding eigenvector. 
After substitution into Eq. 3, the generalized eigenvalue 
oroblem becomes, 

Kd = w°M@ (ai. (5) 
Written in ere form for several eigenvalues, Eq. 5 
becomes, 


Kb = 2°Mo (Eq. 6) 


ies 








Where $ isthe modal matrix, and 2* = diag (w?) is the spectral 
matrix. 

SEW elemsckcdire analysis, tie wqilobalsseiteness 
matrix is already calculated, then only the global mass 
matrix evaluation is needed. Both generalized and lumped 


mass options are coded. These are defined respectively as: 


NE Oe 
oo (Eq. 7) 
i=l 6 I, ZI, 
and 
NE 
M = 2 Sahivi [| (Eq. 8) 
- i=l 6 6x6 


Where I is the identity matrix, the sub-space iteration 
method of Bathe and Wilson [Ref. 9] is used to solve for a 
specified number of lowest eigenvalues and the associated 
eigenvectors. The method is economically efficient for 
large problems. The mass matrix may be diagonal or banded. 
The method is well suited for re-analysis when small 


changes are made in the design. 


C. ANALYTIC GRADIENTS OF THE CONSTRAINTS 

The necessity to compute gradients of the relevent 
functions in a design optimization process arise from the 
fact that efficient mathematical programming algorithms 
require information on derivatives. Furthermore, approxi- 
mate methods based on a Taylor Series expansion of functions, 


i 





requires determination of derivatives. Based on the static 


and dynamic analyses, the gradients of forces, displacements, 


and frequencies with respect to the reciprocal of the cross- 


sectional areas and the coordinate variables, 


Ps 


also 


Gradient of Member Stresses with Reseece seo) ene 
Reciprocal of Area Variables 


Stress in a member is defined as: 





a4 
ai ae (Eq. 9) 
ale 
Rents: 
eel 
a = A, (egw Le) 


The partial derivative is then: 


a0. 











25 1 Ri 9 
xe = ox | ae a me oe ECE 
~e ~e a iS 
The first term of the right hand side is zero, then 
S)0/o » ae 
=x = = os (u.) (Ege ee} 
ac din, 
This can be written in exrlicit form as 
T 
00: BD is 
~ = = = ax (a; ) (Eq. 13) 
ES tg 


Ze 


Gradient of Nodal Displacements with Respect to the 
Reciprocal of Area Variable 


Consider the equation 


Ku =f (Eq. 14) 


~ 
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are formulated. 





The derivative of [K u] with respect to some variable X., 
1s 


3 0 
Zz [Ku] = 3x7 (E) (Eq. 15) 
a au 


In this case Xo = 1/A, and the loads are constant and in- 
dependent of the areas (weight of the truss elements are 


ignored), then 


~ 2 ~Q 
@s 
K Bae = Sey (E 17) 
Sox, ae ae 
~ ~~ 2 
Finally, then 
=~ (4) = ee seo (Kyu (Bea eeks) 
aa - ~ 2 


where K+ 1s the inverse of K. It 1s necessary to compute 
the partial derivative of the global stiffness matrix K 
defined as: 

E. 


NE py 
ee te Te 
i=l 2 





where D; is the matrix of direction cosines defined earlier. 
Therefore, the partial derivative with respect to the 


reciprocal of the area variables is defined as: 


NE EA; 
(K) = > bo [D. ] (Eq. 20) 
i=l 





Be 
L 


Qo” 
2d Qo 
2 
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where 95 1s the kronecker delta defined by o. =lifi= ez 


2 
and Oso = 0 if 1 # 2. In practice the full matrix, SX 
2 
1s not actually stored. Instead the product, se (K)u 1s 
oo 


created directly. Details of efficient gradient computations 


a 





(K) 


as 


are given in [Ref. 10]. 


3. Gradient of Beeguem@encs Wo cee ope CEE the Recep nocal 
of Area Variable 


Consider the eigen-problem defined by, 
ce w oi] o. = 0 (Eq. 21) 
Taking the derivative with respect to the variable X [Ref. 11] 
gives 


ae ie 


- (Eq. 22) 
9 


wien 
1 
= 
7 
= 
oQ- 
_— 
il 
© 


then 


<i {{K - w2M]} 9, + [K - wiM] eo, = 0 (Eq. 23) 


IX 


Pre-multiplying through by b and applying the condition of 


symmetry of the matrix Es ~ aM] gives 





3K OM dw. 
ee eee arent) a Boa ox, ea xe nM Os 
(Eq. 24) 


The left hand side of the equation is zero because of Eq. 2l, 


and thus aK , «OM 
2 gh 7 iia? 
ae ale ee (Eq. 25) 
ee it 
ae oes 
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The partial derivative of the generalized mass matrix with 


respect to the reciprocal of the area variable is defined as: 


3 


es 


(EG 26) 


Qo 


HN 


t 


NE : 
fe OA a 
i=l 


The partial derivative of the lumped mass matrix with respect 
to the reciprocal of the area variable is formulated following 
the same procedure as the generalized mass matrix. 


4. Gradient of Stresses with Respect to the Joint 
Coordinate Variables 


Srieceouseress in a member is 





Bivi 
a 
the gradient is calculated as 
a0. au 
ao: 3 ll At aL 
wo FE we (CS K.) Ut eK. OS Ciclo 
IX, IX, A.=i i A. a: 3X 


In this case, the stiffness matrix is a function of the 
coordinate variables so the first term on the right hand 
Side of Eq. 28 is not zero. 


5. Gradient of Displacements with Respect to the Joint 
Coordinate Variables 


Consider the following equation: 


u=F (Eq. 29) 


~ 


vn 


The gradient of u with respect to xX deSe- 


eG) = -K7t 9 


K (K)u (Eq. 30) 
ax, S 35, °3 





The partial derivative of the global stiffness matrix K now 


1s defined as: 


3 — (AGE, 28; 3AQEy : 
se (Ky =) (RE RE D (-1) 3 (ax) [K,) 
IX = L2 Xx 2 ae a 
1=1 a 
(iar. 3) 
where (-1)/7 and j} = 1 or 2 defines the sign of the gradient 


at that particular node. Note that the terms in the summa- 

tion are only evaluated for members connected to the joint 

defined by the variable Xo. 

6. Gradient of the Natural Frequencies with Respect to 
the Joint Coordinate Variables 


The value of the derivative of the eigenvalue w? 


found from 





epee = 
[K weM]o; 0 (cys 2 ) 
iS 
KOM 

ee | Ske > aces 

os af Mo y (E 33) 
ax, T 4° 

=i vy OF 


The partial derivative of the generalized mass matrix with 


respect to the coordinate variables are defined as: 


NE 
aM 7 Omnpcrem 2 1. 
ee, — a > (anes ee ei (EG 34) 
IX, 6L ela 
i=l a 
where (-1)/ and joe lOone 2 decineomenewsion or the gradient 


at the particular node. The partial derivative of the lumped 


mass matrix follows the same procedure. 
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D. APPROXIMATION CONCEPTS 

At this stage the analysis tools necessary for the design 
process under static and dynamic conditions may be regarded 
as available. It has been pointed out that the application 
of mathematical programming methods for structural design is 
the most widely used because of its great generality and its 
Simple formulation. However, it is required that a large 
number of structural analyses and sensitivity analyses be 
performed. This has motivated the idea of formulating simple 
and explicit approximations for the most relevant response 
quantities. 

These approximations can only be expected to be of 
acceptable quality in some finite region of the design space 
Surrounding the point about which the approximations were 
constructed. The total number of analyses required to find 
an optimum design uSing approximation concepts 1S significant- 
ly less than the number previously required. 

Some sources of Simplification can be considered. First, 
the dimension of the design space can be reduced if a proper 
subspace can be identified. Second, linking of design 
variables which is imposed because of symmetry or practical 
considerations also reduces the dimension of the design 
space. 

The objective function for the design of trusses is a 


Relemively simple explicit function of the design variables. 
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On the other hand nonlinear constraint functions are very 
complicated. From a computational standpoint, a small 
portion of the constraints play an active role in the opti- 
mization process; therefore, deletion of non-critical 
constraints avoids effort of evaluation of irrelevant 
constraints. 

A method to deal with complicated constraint which is 
very effective in reducing computational effort is to deduce 
Simple and explicit expressions for the constraints. 
Linearization is directly and efficiently accomplished by 
Taylor series expansions. The order of the expansion 
selected is decided based on the degree of nonlinearity and 
the approximation required; and a trade-off must be made 
between the computational effort required for the highest 
order derivatives versus improvement of approximation. 

Application of a first order Taylor series expansion on 
stress constraints has been found to be sufficient. This 
was not the case for the natural frequency constraints, 
which may still be numerically unstable during the optimi- 
zation. This suggests further expansion up to a second 
order. However, the computational effort needed to do this 
usually exceeds the benefits. Therefore, a first order 
expansion with move limits is used. 

First order Taylor series expansion of a function, W, 
Meeiemcerdole sx, about a point, XO, 1S written as 


W(p)= W(x) = W(Xo) + (X-Xo) V W(Xo) 
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The philosophy underlying the use of linear approxima- 
tions 1S not to transform the problem into a sequence of 
linear programs, but to replace the constraints by simple 
aad explicit approximation functions. In this study the 
Taylor series expansion was used for all constraints with 
respect to the reciprocal of the member sizing variables, A. 
The objective function, in reciprocal space, 1S now non- 


linear but Still explicit and easily evaluated. 


Peers scr iIVE FUNCTION 

The function whose least value is sought in the optimi- 
zation procedure is defined as the total weight of the truss 
which is given by 


NE 


w=) O.A.L. (Eq. 35) 
1=1 


Where NE is the total number of members, and O; is the 
material density, both are prescribed constants. A; and 


L; are the area and length of the ith member respectively. 


ree CONSTRAINTS 

The restrictions to be satisfied in order for the design 
to be acceptable are formulated explicitly; behavioral and 
side constraints are defined accordingly. 

iS GCSiS 


; EG. 36 
[ei £95 eS Obi (Eq 
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where Ong and Oo,; are the allowable compressive and tensile 
pereosoes respectively for member i. 
Zemeeuler Buckling 


The Euler buckling stress is formulated as 





_ CA.E, 
where c 1S a prescribed constant. 
3. Frequency 
fc ans 
Wag Wag Wr GEG. 36) 


= ~ ¢ 
where Wr and W,, are the lower and upper bounds respectively 
on the first natural frequency of the system. 

4. Limits on Areas 


< < 
“min; a Ay ~ “max. ge 22) 


where, Bee poand Bare are the minimum and maximum allowable 
Z i 


cross-sectional areas of the ith member, and are taken to 
be the same for all members. When symmetry of deSign is to 
be preserved, linking of the variables is required. This 
is defined as 


A, = A; (Eq. 40) 


where K and i are symmetric members. 
5. Displacement 
~ + 
S S ~ 41 
U Su, S uy (Eq ) 


where u,, and u. are the lower and upper bounds respectively 


K K 


Smereiit displacements Uy 
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6. Limits on Coordinate Variables 


- + 
Xi is x < Ky (Eg. 42) 


where again oe and xe are the lower and upper bounds 
respectively on the Kth coordinate variable, and if Symmetry 
must be preserved, linking of thewariables is required. 
This is obtained by 

X, = a, + b, Xe (formas) 


where ay and Dy are constants and X., 1s the coordinate 


K 
Variable. 


G. GENERAL FORMULATION 
The inequality constraints are of the form, 
G5 (X) 0, J = eee NEC (eae) 


The constraints are normalized as follows: 


DEEeCSSes act eG 
Pt 
een (Eq. 45) 
6 

Euler Buckling eo (Eq. 46) 
Ob 


Note that because o and o), depend on the member area A,, 
this constraint is treated as a nonlinear function. The 
Taylor series expansion is performed on o and the value 


Or 3}, is continually updated. 
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LA 


Displacements os = 0 
6=— 
a 0 (Eq. 47) 
6+ 
rrequency “ 5156 
W) 
The side constraints are 
< 
Ay = A, - 0 (ge 49) 
ee 3 
Ay Ay 0 
xr - x, $0 (Eq. 50) 
K K 
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wo LNT RODUCTION 

The main goal of structural engineering is to design 
structural systems that efficiently perform specified 
purposes. Selection of a specific algorithm is required 
and this algorithm must minimize the number of times the 
Structure has to be analyzed and the amount of specific 
gradient information required. Finally, the algorithm 
should provide reasonable assurance that it will attain an 
optimum or near-optimum design. 

The next two sections are a brief explanation of the 


algorithm used for this work. 


B. GENERAL FORMULATION 
The general constrained minimization design problem is 


defined as 


minimize 

W(X) (Eq. 1) 
Subject tO 

G(X) <0 een (Eq. 2) 


where, W(X) is the objective function. Functions G(X) are 
the set of inequality constraints. The vector of design 
variables X includes member sizing variables a and geo- 


metric design variables Kae 
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C. OPTIMUM GEOMETRY DESIGN 

The procedure used here was to treat the geometric design 
parameters as independent design variables. The member sizing 
parameters are handled as dependent variables which are 
determined as a sub-problem. 

Beginning with an initial geometric design vector Xr the 
design proceeds iteratively using the following relationship: 


aa a x 
yor = xo + a Sq 
g g 9 g 


teas. 23) 
where q is the iteration number and Sy is the search direction 
to be determined. a is the scalar parameter determining the 
distance of travel in the design space. 

For each proposed geometric vector, Kor the structure 
was optimized with respect to the member sizing variables, 


Xa! 


by the sub-optimization problem defined in Section D. 
Assume that for the initial geometry the structure has 
been optimized with respect to the cross-sectional areas, 

and that, from this subproblem there are 2% active constraints 
Serene LOrm: 

G(X) ed k = Wee. £ (Eq. 4) 
where, G(X) is defined as active if its value 1s close to 
zero. 

Now, it is necessary to find the search direction, So) 
so that by moving in this direction in the coordinate 


design space, the objective function is minimized. This 


direction is found by solving the following subproblem. 
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Minimize 


VF(X).S (Eq. 5) 

Subject to: 
G(X) wo = Sane: ogee GHGa 6) 
S.s<1 (Eq. 7) 


At this point, S, provides a search direction in the 
combined [X, + x] Space which is projected onto the active 
constraints. Only the S part of S is used. 

Once, Sor is known it is substituted back into Eq. 3 
and a one-dimensional search on a 1s performed to update 
the vector, Xo The Optimum-Geometry design problem can be 
summarized in the following algorithm. 

1. Given an initial coordinate design vector, Xo and 


area design vector, A, - SPCC LEv, a» and Be 


max ax’ (eae! 


maximum change in the coordinate at each iteration, and 


Pp total number of iterations.) 


max’ 


2. Solve the fixed geometry problem and calculate the 
minimum weight (W*) for the current geometry. 
3. Determine the set of active constraints. 


4. Determine the search direction vector, S. 


* 


5. Find the move parameter chy On ua De change of some 


ax 
coordinate variable, or such that some coordinate constraints 
becomes active. 

6. Solve the fixed geometry problem for coordinates. 


2 - k= = 
Ro = Ag + A156 and calculate the optimum weight, W(X). 
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7. Find the parameter ae which minimizes the weight 
using 2-point quadratic interpolation. If aL >a, do to 
step 9. 

8. Solve the fixed geometry design problem for coordi- 
nates xa = XS + sé and calculate the optimal weight. 


9. Check convergence; if satisfied, terminate; otherwise, 


set q = qt+tl, update X and return to step 3. 


Pee ve GEOMETRY DESIGN 
As stated earlier, a sub-optimization problem has to be 
solved for a proposed geometry design vector, X. The 
structure 1S now optimized based on the cross-sectional 
area subspace. This is defined in general as: 
Minimize, 
W(X) ( Bepeeeo) 
Subject to, 
G(X) aa) = le eee (Eq. 9) 
The design proceeds iteratively. Given an initial vector 
of design variables Xo find the X vector at the (qt+tl) ith 


iteration defined as: 


om Jt m= 
xatl - xd 4 x" 89 (Eq. 10) 
m m m m 


* e = om ° 
where, Ca is a scaler multiplier, and, San’ 1S a vector move 
direction in the design space. 
Now the problem becomes one of finding the direction, 


S, and the move parameter, a*. Zoutendijk [Ref. 12] shows 


34 








that the direction may be found by solving the following 
problem. 
Maximize 
8 (qe ay) 


Subject to 


Tse) aS 4 0 (Eq. 12) 
eee js 5 + 6.8 S = er , 

ai =! 4° 0 5 i NAC (ache ins) 
Se Suse) (Eq. 14) 


where the scalars . are named as "push-off" factors. 
If some of the constraints are violated, this algorithm 
is modified in order to find a feasible design [Ref. 13]. 
For the fixed geometry sub-problem, the reciprocals of 
the member sizes are used as design variables and approxima- 
tion techniques are employed [Ref. 5]. The optimization 


program CONMIN [Ref. 14] is used to solve this sub-problem. 


35 


an 2 





IV. NUMERICAL EXAMPLES 


A. INTRODUCTION 

Design of planar trusses and space towers are presented 
here and the corresponding numerical results are summarized 
to demonstrate the purpose of this study. 

The examples begin with an 18-bar truss. In this case, 
it 1s shown how the optimum geometry is dependent on the 
constraints imposed. A single load condition was considered 
and the cross-sectional areas are linked. 

Next, a 47-bar planar tower is designed to support a set 
of load conditions given in Table VI. The design is subject 
to constraints on the member stresses, Euler buckling, dis- 
placement and first natural frequency. Linking was imposed 
for symmetry in both cross-sectional area and coordinate 
variables. 

Finally, a 25-bar space tower was designed. Constraints 
on stresses, Euler buckling, displacement, and frequencies 
were imposed. The truss was required to support two 
different load conditions. Symmetry for the member areas 


and coordinates was established. 


B. CASE 1: 18-BAR PLANAR TRUSS 
A cantilever truss, as shown in Fig. A.1l, has been used 


previously as an example for the design of trusses of a 
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specified geometry [Ref. 7]. This structure was designed 
for optimum geometry subject to a Single set of load condi- 
tions given in Table I. The allowable stresses are specli- 
fied as 


-20,000 < Oo; <0, ,000m ps1 


Young's modulus is taken as 10’ pS1 and the material density 
09 = 0.1 1lb.fcu. in. The allowable stress at which Euler 


mack ling Occurs is 


fee 


The independent coordinate variables were taken as X3, 
Y¥3, X5, Y5, X7, U7, X9, Y9. The member areas were linked 
in the following groups: Al1=A4=A8=A12=A16; A2=A6=A10=A14=A18; 
A3=A7=A11=A15; AS=A9=A13=Al7. There are a total of eight 
independent coordinate variables and four independent area 
variables. 

1. Case la 

The structure was designed subject to stress constraints 

only. The resulting geometry is shown in Fig. la, and the 
design information is given in Table II. Weight versus 
iteration history is plotted in Fig. A.2. The number of 
analyses for this design is 59, and the execution time 2.53 


seconds on an IBM 3033 computer. 
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2. Case lb 
stress and Euler buckling constraints were imposed 
for this case. Design information is given in Table III. 
Figure lb shows the final geometry layout and the weight 
versus iteration histcry is shown in Fig. A.3. The number 
of analyses for this design is 78 and the execution time 
was 3.94 seconds. 
3. Case lc 
This design is based on stress, Euler buckling and 
displacement constraints. The latter was applied at node 
number 1 in the Y-direction. The results are shown in 
Table IV. Figure lc presents the final layout and the 
weight versus iteration is plotted in Fig. A.4. The number 
of analyses is 91 and the execution time was 3.94 seconds. 
4. Case ld 
This final case includes all the constraints mentioned 
before plus a constraint that the first natural frequency 
be greater than or equal to 3 Hz. A non-structural mass 
of fixed value W = 1,000 lbs was placed at node l. Steady 
convergence is achieved and results are summarized in Table V. 
The final geometry and weight iteration history are shown in 
Figures ld and A.5, respectively. The number of analyses is 
96 and the execution time was 4.73 seconds. 
When the mass is removed from the structure, the 
design fails to converge even when move limits are imposed. 
First, it is known by definition that the frequency is 
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proportional to the stiffness and is inversely proportional 
to the mass. Second, both the K and M matrices are functions 
of area and coordinate variables. Therefore, as K increases, 
M also increases, and the frequency is kept closed to the 


initial range. 


C. CASE 2: 47-BAR PLANAR TOWER 
The initial layout of the tower is shown in Fig. A.6. 
stress constraints were imposed as well as constraints on 


Euler buckling, displacement, and first natural frequency: 


-15,000 s Oo; S207 080n psi 1=1,47 
and 
10.17EA, 
Oo. 2 0,. = - —— —— psi 
: ba gL 
a: 
ee 12 hz. 


respectively. The members are assumed tubular with D/t=10. 
moung S modulus of 3.x Toe psi, and material density 

R= .3 lbs/fcu. in. Minimum allowable area of We in. was 
imposed. Symmetry about the y-axis was desired during 
optimization so linking of variables is necessary. Joints 
15, 16, 17, 22 are kept fixed, and joints 1 and 2 allowed 
to move in the x-axis direction (Y=0). This gave a total 
of 27 area design variables and 17 coordinate variables. 
Nonstructural masses were attached at nodes 19 and 20, of 


W=500 lbs each. The results are tabulated in Tables VII 


through XII. The final geometry for the case where all 


a2 





constraints were imposed is shown in Pig. ance the 
Bectaeton history is shown in Figs. A.8 through A.10. When 
all constraints are imposed, the design required 220 analysis 


and the execution time was 160.7 seconds. 


D. CASE 3: 25-BAR SPACE TOWER 

The 25-bar space tower shown in Fig. A.1l was designed 
to support two independent load conditions given in Table 
XIII. The allowable stresses were specified as 


-40,000 < a; S$ 40,000 psi 1=1,25 


Young's modulus was selected as 10/ pSi and the material 
density o= .1 1b./cu. in. The members are assumed to be 
SimoulaG with a nominal diameter to thickness ratio of 
D/t = 100., so that the stress at which Euler buckling occurs 
1s 

LO at EA; 


1. = - ——— psi. 
ioe gL 


Symmetry with respect to both x-z plane and the y-z plane 

was imposed, so linking of variables were made as follows: 

the member areasS were grouped in the following sequence: 

Al, A2=A3=A4=A5, A6=A7=A8=A9, A1LO=All1l, A12=A13, A14=A15=A16=A17, 
A18=A19=A20=A21, and A22=A23=A24=A25. For the coordinates 
wee A428 8G, 18 were Considered aS variables. The joints 

1 and 2 were held fixed and joints7 through 10 were required 

to lie in the x-z plane. Nonstructural masses of W=500 lbs. 
were attached at nodes 1 and 2, respectively. The first 
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natural frequency was limited to a value \ 2 16 Hz. 
Results of this example are shown in Table XIV and the 


heermation Nistory is shown in Fig. A.12. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The final layout dependency on constraints has been 
presented for the elastic design of trusses for optimal 
geometry. The truss may be planar or three-dimensional and 
may be indeterminate. Multiple load conditions were con- 
Sidered. The design procedure was separated as: analysis, 
design for fixed geometry, design for optimum geometry. 

The displacement method for static analysis and the 
subspace iteration method for dynamic analysis were applied. 
The sequential optimization based on two design sub- 

Spaces present substantial advantages in the reduced 
number of analyses and allow the designer to keep control 
of the optimization process. 

Several examples were considered. In every case, signifi- 
cant weight reduction was efficiently achieved. Also, the 
geometries obtained appear quite acceptable from an 
aesthetic as well as structural point of view. 

The graphs of Weight vs. Iteration number show that an 


acceptable design can be achieved in few iterations. 
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B. RECOMMENDATIONS 

The following recommendations may be of theoretical 
and practical value. 

1. Weight due to the structure itself and other design 
dependent external forces can easily be taken into account 
and should be considered in future studies. 

2. Application should be made to reasonable sized 
structures such as offshore towers and long span roof 
trusses. 

3. The principles and the procedure described herein 
can also be used for optimal design of frames as well as 


other structures. 
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APPENDIX A 





Figure 1. 18-Bar Planar Truss. Initial Geometry. 


Wi= 4,780. (lbs) 


We 4,524. (lbs) 


Figure la. Stress. 


Wi= 6,430. (lbs) 
Whe 5,724. (bs) 


Figure lb. Stress, Euler Buckling. 


Wis 7,413. (Tbs) 


Wfz 6,465. (lbs) 


Figure le. Stress, Euler Buckling, Displacement. 


Wi= 13,686. (lbs) 
Wf= 11,528. (Ibs) 


A 
Figure ld. Stress, Euler Buckling, Displacement, and 


Frequency. 


Figure A.1. 18-Bar Planar Truss. Stress, Euler Buckling, 
Displacement, and Frequency Constraint. 
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Figure A.6. 47-Bar Planar Tower. 
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TABL 


18-Bar EFlanar Truss. 


LOAL CCNDITION 


Fx 
Q. 
Q. 
is 
Q. 
0. 


Avestalee 


oO OH £ hd 


ee 


boadsS and Constants. 


TREE ChE Cote ANTS 


Younea* oetodu lus E 
Allewatle Stress Tr 
Density i 


Buckling Goecrricient 


C= 


(lbs.) 
Fy 
eZ OO 
ZOO Oe 
OOOO. 
ee ONG. 
eg OU 
Q. 1EC8 isl. 
OR eed Ss) ecole 
OSes / Gus h 
~4E08 Pol. 
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(.69E08 KN/m ) 
(.138E06 KN/n ) 
(.276E04 Kg/n ) 
(.276E0S KN/m ) 





TABLE IT 


[S-Bearcwbtanar Truss. Pesign PerOrmation. Stress: 
ARE A ™® ({sq.in.) 

Member Initial Fano 
A1l=A4=A8=A12=A416 10.0 11505 
AZ2=A6=A10=A14=A18 1520 15 07 
A2=A7=A11=A15 5.0 4.54 
A4=A9=A13=A17 7.2.07 S35 

CCORCDINATES (in.) 
Jean Initial Final 
x Yy x Y 

3 1000 0.0 991.17 19.686 

5 750 0.0 745.88 152 12 

7 500 G0 494.61 S07 Sac 

9 250 0.0 249.54 23.609 


Cetinum Weight for Initial Geometry 


Cetimumn Weight for Pin 


al Gecmetry 


4,780.5 (lks.) 
4,524.7 (lbs.) 


* Areas are the optimum values fer the initial and final gecmetry. 


a 





Tabi ELL 


18-Bar Planar TIruss. Design Informaticn. 


Buiextcsnug. 


AOR EB AS®.(sg.in.) 


Member 


A1=A4=A8=A12=A16 
A2=A6=A10=A14=A 18 
A3=A7=A11=A15 
A4=A9=A13=A17 


Initial 


10.00 

21.65 

2a 
7.071 


CCORDINATES (in.) 


JEint 
x 
3 1000 
5 750 
7 500 
9 250 


Tonat 1a lL 


Cetimum Weight fer Initial Gecmetry 


Crtimum Weight for Final Geometry 


* Areas are tke optimum values for the initial and 
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Stress, Euder 


Final 

11.34 

19.28 

10.97 

5.306 
Final 
4 7 

994.57 162..3 1 
747.36 10229": 
482.90 ize Jo 
227 172105 


6,430.7 (lks.) 
5,713.C (lbs.) 


final gecmetry. 





18—-~Bar Elanar Iruss. 


TAREE IV 


Design Informaticn. 


Buckling, Displacement. 


Member 


A1=A4=A8=A12=A16 
AZ=A6=A10=A14=A 18 
A3=A7=A11=A15 
A4Y=A9=A13=A17 


gCrne 


Oo ws Wm ty 


APR YE UA * o(Sq- 20.) 


Initial 


14. 18 
21.66 
TZ, 5 

10522 


COOREINATES (in.) 


1000 
750 
500 
250 


Patras 


Cpetimum Weight fer Initial Gecmetry 
Cetimum Weight for Final Gecnretry 


Stress, Euler 


Final 

16.27 

20.06 

171.18 

7.863 
Final 
x y 

962.19 109.93 
703.95 55.3309 
452.7 215307 


= 7,413.1 (lbs.) 


6,466.9 (lbs.) 


* Areas are the optimum values for the initia] and final gcemetry. 


a 





TABLE V 


18-Bar Planar TIruss. Design Informaticn. 


SBcLess, Euler 


Buckling, Displacement, Frequency. 


Member 


A 1=A4=A8=A12=A16 
A2=A6=A10=A14=A 18 
A3=A7=A171=A15 
A4=A9=A13=A17 


AR 


EAS 


(sgosn.) 


Ppatial 


21.68 
40.58 
12.5 
Bo \euo 


CCCRDINATES (in.) 


Jcint 


Ww 


cn 


wo -~ J 


Cetinun Weight fer 
Crtimum Weight for 


* Areas are the optimum values 


1000 
750 
500 
250 


Initial 


Initial Gecmetry 


Final Gecmetry 


Final 


25.66 
32629 
13.03 
16.09 


Final 


G7i2a52 
717.08 
46€.75 
200.98 


88.0 

27.84 
Seo | 
m9, 29 


1o7090.0 (LES. } 
tice. c (lbs-) 


for the initial and final gecmetry. 
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TABLE VI 


4“7-Barc Planar Tower. 


LOAD CCNDITICN 1 
TeCINT PX 
le, 6,C00. 
D0 O. 
LOAD CONDITION 2 
Joint FX 
17 0. 
Z2 6,000. 
LOAD CCNDITICN 3 
Joint FX 
17 6, C00. 
22 6,000. 
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Road. Cond2 <2 Oonse 


(lbs.) 
Fy 


-14,000. 
QO. 


(lbs.) 


FY 


-14,000. 


(ik. ) 


PY 


-14,000. 
-14,000. 





TABEE ov EL 


iv Bac Planar Lowes. Azea. Stress, Euler Buckling. 


Menrker inetiad Pinal 
Area * (sq.in.) Area * (Sqg.in.) 
3 3.764 2.727 
4 3.315 2.468 
5 05787 Oe 2d 
v] 0.864 Oe2i3 
8 0.856 0.938 
10 1.754 1.076 
12 2.087 1.691 
14 1.188 0.655 
is ozo 1.058 
18 2.007 1.412 
20 0.648 02263 
22 0.843 Ome A 
24 1.700 1.060 
26 1.700 T0s2 
27 1354 0.820 
28 0.847 0.302 
30 3.609 2.766 
31 list 5 0.657 
33 0.638 0.207 
35 2.842 Ze cod 
36 0.676 0.266 
38 1.596 1.408 
4Q 3.686 3.429 
41 1.526 0.991 
43 0.677 0.170 
45S 4.486 2650 
46 dies 2 1.005 


* Crtimum areas for initial and final geometry. 
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(°sqt) t26°2 = Kx,yawmoayH Teurtg AoZJ rybtamg_ anatydo 


(esqt) €06’°2 = Axr29mo0aH TeTITUT TOF yybrax wnatydo 


0°009 o°O0St 0°909 "OSt Gc 
LtL°eL9 £0°66 0°009 "06 A 
hB°SZ9 GZ°6t 0°009 "OE 02 
0°OLS 0°06 0°OLS "06 9t 
ZE°6LS 0°09 0°OHS “OE he 
6S°08h 60°6S 0°08h GS aot 
EL Ozh 0°09 0°0Zh "OE | OL 

B9°LSE BB°L9I Owe "09 
99°6EZ 88°RL - oon “09 9 
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TABLE X 
47-Bar Planaz Tower. Aréas. Stress, Fuler Buckling, 


Displacement. 


Memker Initial Final 


Area *® (sg.in.) Area * (sq.in.} 


3 15.28 4.478 
4 ans 3.671 
5 3.066 1.037 
7 2.691 0.561 
g ASE 1.156 
10 3.706 1.170 
12 Sez 2esa6 
14 4.048 0.€69 
15 S000 1.385 
18 4.780 1.169 
20 1.650 0.562 
22 26233 1ezigs 
24 4.665 1.526 
26 4.816 15544 
27 3.693 1.036 
28 2.923 0.566 
30 17.91 San 
31 7.091 0.949 
33 2.566 0.558 
SG igs Ses 
36 6.640 tat? 
38 3.385 1.076 
40 aaah Tao 
41 6.792 ie 70 
43 2-708 Gmecs 
4S 37.11 9.878 
46 7.363 Ze 


* Optimum ares for initial ard final geometry. 





eA is ek 
47-Bar Flanar Tower. Areas. Stress, Euler Buckling, 


Cisplacement,and Frequency. 


Meurer Initial Pinal 
Area * (sq.irz.) Area ® (Sq.in.) 
3 16.59 Ue Z5 
4 15.30 4.6065 
5 V2 | 0.938 
a eods 0.348 
8 to3o6 1.081 
10 Sh lise i039 
12 acoo3 3.025 
14 5.988 1.088 
15 &.468 1.6951 
18 8.593 V.656 
20 1.416 02645 
22 1.836 1.724 
24 8.470 2.480 
26 8.470 2.186 
27 5.999 tea 
Ze 1.949 0.452 
30 14,54 u2 233 
31 7.096 0.855 
33 1.709 O37 
35 9.710 4.579 
36 1.414 0.368 
38 1.385 1.439 
40 10.43 4.396 
4 1.484 1.126 
43 1.465 0.270 
45 Wie 4.365 
46 1.482 1.664 


* Crtinum areas for initial and final geometry. 
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TACC Ee LY 


25-Bar Space Tower. Design 


COORDINATES (in.) 


Intormetion. 


Joint Taietial Strs-&ck. Displacement Prequency 
x a7 <3 21.487 Jae a2 64.9713 
4 y 37.5 4&.271 602.83 47.271 
z 100. 100.27 114.26 128.77 
; x 100. 22.146 47.383 107556 
y 100. 96.353 140.53 110.077 
AR EASs@en(SG. on.) 
BENBER STRESS~BOCKL CISELACEMENT PREQUENCY 
INITIAL FINAL INITIAL FINAL INITIAL FINAL 
1 0.9E=+3 0.013 0.668 0.029 0.548 0.070 
2 0.782 0.414 2.032 0.450 Zeo34 0.73% 
6 0.754 0.842 3.100 1.219 2635 lez 
10 0. 1E=3 Dos 0.533 0.015 0.547 0.054 
12 0.130 a 1a 0.549 0.118 C.496 0.206 
14 0.558 0.121 0.668 0.084 Oreo 34 0.350 
18 0.982 0.739 los 2 0.751 1.449 0.720 
22 0.801 0.554 £2003 0.668 3.363 1.130 
Stress-Buckl. Cisplacement Frequency 
initial Pinal Teitral Franal Toitsat Final 
R@Ight (ib<.) 229.52 128.55 2605.2 169.9 SB) 473 261.52 
Iteration 12 44 9 
Tise (sec.) 5.38 37.14 21.09 


* Optimua areas for initial and final geonetry. 
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